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Abstract 

The ability to tailor the porosity and optical properties of columnar thin films (CTFs) renders them 
promising platforms for optical sensing. In particular, surface-plasmon-polariton (SPP) waves, guided by 
the planar interface of an infiltrated CTF and a thin layer of metal, may be harnessed to detect substances 
that penetrate the void regions in between the columns of a CTF. This scenario was investigated theoretically 
using a higher-order homogenization technique, based on an extended version of the second-order strong- 
permittivity-fluctuation theory, which takes into account the size of the component particles which make up 
the infiltrated CTF and the statistical distribution of these particles. Our numerical studies revealed that as 
the size of the component particles increases and as the correlation length that characterizes their distribution 
increases: (i) the phase speed of the SPP wave decreases and the SPP wave's attenuation increases; (ii) the 
SPP wave's penetration into the CTF decreases; (iii) the angle of incidence required to excite the SPP wave 
in a modified Kretschmann configuration increases; (iv) the sharpness of the SPP trough in the graph of 
reflectance versus angle of incidence increases; and (v) the sensitivity to changes in refractive index of the 
infiltrating fluid decreases. 

Keywords: Surface plasmon polariton, columnar thin film, strong-permittivity-fluctuation theory, Brugge- 
man formalism 

1 Introduction 

A columnar thin film (CTF) consists of a parallel array of columnar nanowircs, oriented at angle ^ to a 
planar substrate [1]. The columns may be grown on the substrate by means of physical vapour deposition. 
Through judicious control of the vapour deposition process, both the macroscopic optical properties and the 
porosity of the CTF can be tailored to order [2] . Consequently, CTFs are promising candidates as platforms 
for optical sensing applications, wherein it is envisaged that the species to be sensed is contained within a 
fluid which penetrates the void regions in between the CTF's columns. 

A recent theoretical investigation highlighted the fact that surface-plasmon-polariton (SPP) waves, guided 
by the planar interface of an infiltrated CTF and a thin layer of metal, may be usefully harnessed for optical 
sensing [3]. In this scenario, the angle of incidence required to excite an SPP wave was found to be acutely 
sensitive to changes in refractive index ni of the fluid which fills the CTF's void regions. Furthermore, the 
phase speed and propagation length of the SPP wave were also sensitive to ni. 
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A crucial step in modelling CTFs — as well as more general sculptured thin films — as platforms 
for such optical sensing is the estimation of the effective permittivity dyadic of the infiltrated CTF [4]. 
This may be achieved in a two-step homogenization process: First, an inverse homogenization procedure 
can be implemented to estimate certain nanoscale parameters (namely, refractive index, shape parameter 
and volume fraction) for the uninfiltrated CTF, using experimentally-determined values of the uninfiltrated 
CTF's permittivity dyadic [5]. Second, a forward homogenization approach can then be used to estimate 
the effective permittivity dyadic of the infiltrated CTF. In recent studies the well-established Bruggeman 
homogenization formalism was employed for this purpose [3, 6, 7]. However, there are certain limitations 
with the Bruggeman formalism. It does not take into account the nonzero size of particles which make 
up the component materials nor does it take into account the statistical distribution of these component 
particles. But both the size and distribution of the component particles can have a significant bearing 
upon the estimates provided by homogenization formalisms. For example, significant higher-order effects 
associated with the nonzero size of component particles were recently reported for a homogenization study 
based on an infiltrated chiral sculptured thin film [8] . 

In the following we investigate CTFs as platforms for SPP optical sensing using a higher-order approach 
to homogenization, based on the strong-pcrmittivity-fluctuation theory (SPFT) [9]. A second-order imple- 
mentation is used wherein the statistical distributions of the component particles are characterized in terms 
of a two-point correlation function and its correlation length L [10, 11, 12]. Also, an extended implemen- 
tation is used wherein the nonzero size of component particles is accommodated via a size parameter 77 
which originates from the corresponding depolarization dyadics [13]. Two SPP scenarios are investigated: 
(i) a canonical scenario wherein the SPP wave is guided by the planar interface of a CTF half-space and a 
metal half-space (Sec. 3); and (ii) a realistic scenario for SPP wave excitation based on a modification of 
the Kretschmann configuration (Sec. 4). Our particular goal is to elucidate the effects of the higher-order 
homogenization parameters L and 77 on SPP wave excitation. As the theory which underpins our study is 
described comprehensively elsewhere, we provide only a brief outline of the theory (along with appropriate 
citations to the recent literature) before presenting our numerical results. 

As regards notation, vectors are written in boldface while 3x3 dyadics arc underlined twice. The unit 
vectors u x , u y and u z are aligned with the Cartesian x, y and z axes. The real and imaginary parts of 
complex quantities are provided by the operators Re [■] and Im[-]; and i = \/— 1. An cxp(— iujt) time- 
dependence is implicitly assumed, with uj being the angular frequency and t being time. The free-space 
wavenumber is denoted as fco- 

2 Homogenization 

2.1 Macroscopic perspective 

Macroscopically, a CTF may be regarded as a homogeneous biaxial dielectric material, characterized by a 
relative permittivity dyadic of the form [1, 2] 



Herein v = 1 denotes an uninfiltrated CTF and v = 2 an infiltrated CTF, while the normal, tangential and 
binormal basis vectors are given by 
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with the column inclination angle x G (0, tt/2] . The permittivity parameters for certain uninfiltrated CTFs 
have been experimentally determined. For example, 
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(3) 



for a CTF made from patinal® titanium oxide, as determined at a free-space wavelength of 633 nm [14]. 
Here Xv is the angle of incidence of the vapour flux during the deposition process used to fabricate the CTF. 
It is related to the column inclination angle x y i a 



tanx = 2.8818 tan^jj 



(4) 



For definiteness, we fix Xv = 30°. From a knowledge of {e a i, £bi, £ci}, the corresponding permittivity param- 
eters for a CTF infiltrated by a fluid of refractive index ng, namely {e a 2, £62, £c2}, can be estimated using a 
inverse/forward homogenization process, as we now describe. 



2.2 Nanostructural perspective 

Within our homogenization approach, the columns of the CTF are envisaged at the nanoscalc as assemblies 
of ellipsoidal particles, all with their major axes aligned with the rotational axis of the columns. Thus, the 
surface of each ellipsoidal particle, relative to its centre, is prescribed by the position vector 

rig-*, (5) 

where f is the radial unit vector which prescribes the surface of the unit sphere. In conformity with eq. (1), 
the shape dyadic U_ has the form 

U = ——= (u„ u„ + 7 r u T u T + 7 h u b u&) . (6) 

A suitably elongated ellipsoidal shape is achieved by choosing the shape parameters 7;, > 1 and j T S> 1. The 
value j T = 15 is selected here because larger values of 7 r do not result in any significant effects for slender 
inclusions [15]. The size parameter 77 provides a measure of the linear dimensions of the ellipsoidal particle. 
In compliance with the homogenization regime, r\ must be smaller than the wavelengths involved, but it 
need not be vanishingly small. In conventional homogenization formalisms, as typified by the Bruggeman 
formalism, the limit rj — > is taken, whereas in extended formalisms nonzero values of 77 are incorporated 
[16]. 

The regions in between the columns of the CTFs (whether infiltrated or not) are taken to have the same 
particulate nanonstructure as the columns themselves. So that there are no spaces in between the ellipsoidal 
component particles, a fractal-like distribution of the component particles is envisaged. Accordingly, the size 
parameter rj should be regarded as an upper bound on the linear dimensions of the component particles [17]. 

In conventional homogenization formalisms, as typified by the Bruggeman formalism, the distributions 
of the component materials are characterized solely in terms of volume fraction: therein, the volume fraction 
of a CTF occupied by columns is / <E (0, 1), while 1 — / is the volume fraction not occupied by columns. 
In addition to volume fraction, the second-order SPFT also characterises the distribution of the component 
materials in terms of a two-point correlation function and its associated correlation length L [9, 12]. The 
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wavelengths involved are presumed to be large in relation to the correlation length. Furthermore, on the 
grounds of physical reasonableness, the inequality L > 77 should be satisfied [13]. 

The deposited material from which the columns of the CTF are composed is assumed to be an isotropic 
dielectric material of refractive index n s . This refractive index is generally somewhat different to the refractive 
index of the bulk material that was evaporated to make the CTF, depending upon the precise details of the 
deposition environment [18, 19]. 

We emphasize that the nanoscale parameters {7^, /, n s } are not readily determined by direct measure- 
ment. Therefore, we rely upon theoretical estimates of these quantities. 

2.3 Inverse homogenization — uninfiltrated CTF 

The first step towards the estimation of the infiltrated CTF's permittivity parameters {e a 2, e&2, e C 2} in- 
volves an application of inverse homogenization. The (noncxtended) Bruggeman formalism is inverted to 
yield estimates of the shape parameter 7&, volume fraction / and refractive index n s , making use of the 
experimentally-determined values of {e a i, e^i, e c i} provided in eqs. (3) for the uninfiltrated CTF. This proce- 
dure is described in detail elsewhere [5]. Parenthetically, certain constitutive parameter regimes are known 
to be problematic for the inverse Bruggeman formalism [20] (and, indeed, also problematic for the forward 
Bruggeman formalism [21]) but these regimes are not the same as those considered here. 

2.4 Forward homogenization — infiltrated CTF 

The second step towards the estimation of the infiltrated CTF's permittivity parameters {e Q 2, £b2, £02} in- 
volves combining the estimated nanoscale parameters {7b, /, n s } with the refractive index of the infiltrating 
fluid ri£ in an application of forward homogenization. In contrast with recent studies, we implement an 
extended version of the second-order SPFT for this purpose. Thereby, the distributional statistics of the 
component materials as well as the nonzero size of the component particles are taken into account. A full 
description of the extended version of the second-order SPFT — appropriate to biaxial homogenized com- 
posite materials, based on oriented ellipsoidal component particles, precisely as is the case here — is available 
elsewhere [13, 17]. Let us note that in the limits L — > and r\ — > 0, the nonextended, zeroth-order SPFT 
emerges which is equivalent to the Bruggeman formalism. Also, higher-order homogenization approaches 
such as those implemented here generally do not scale in a physically-intuitive manner [22] , but the physical 
significance of this aspect is unclear. 

Numerical estimates of e a 2, £b2 and e C 2, as computed using the extended version of the second-order 
SPFT, are plotted against rj/L G (0,1) and kgL € (0,0.2) in Fig. 1. To allow direct comparison with an 
earlier study based on the Bruggeman formalism [3] , the refractive index of the fluid infiltrating the void 
regions of the CTF was fixed at ri£ = 1.5. The real parts of {e a 2, £62, £02} increase modestly, but significantly, 
as both L and r\ increase. More conspicuous is the uniform increase in the imaginary parts of {e Q 2, e&2, £02} 
from the null values that they take in the limits L — > and r\ — > 0. The emergence of these nonzero imaginary 
parts is indicative of losses due to scattering from the macroscopic coherent field [23] . These scattering losses 
increase as the component particles increase in size and as the correlation length increases. Such scattering 
losses are a general feature of higher-order approaches to homogenization, arising whenever the component 
materials are characterized by a (nonzero) length scale [24, 25]. 

3 Canonical boundary- value problem 

Now we turn to the excitation of SPP waves, guided by the planar interface of a CTF half-space and metal 
half-space. This scenario allows us to explore the wavenumbers of the SPP waves and thereby consider 
the phase speeds and propagation lengths of the SPP waves. The theory for this canonical boundary- value 
problem was developed recently by Lakhtakia and Polo [26]. 

Let assume that a metal of relative permittivity e m occupies the half-space z < while a CTF infiltrated 
by a fluid of refractive index ni occupies the half-space z > 0. We consider the p-polarized SPP wave which 
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is guided by the planar interface. The electric field phasor of the SPP wave may be expressed as [26] 



E(r) = ^ 



A m ( u x u z ) exp [ik (ax - iq m z)], 

qm 



iaq c - (e a2 - e b2 ) sinx cosx 
a 2 - (e a2 cos 2 x + e b2 sin 2 x] 



u 2 



z < 



exp [iko (ax + iq c z)], z>0 



(7) 



where A m 
and q m = \fa 2 



and A c are the complex-valued amplitudes, akoU x represents the wave vector of the SPP wave, 
e m . The relative wavenumber a and q c are deduced from the corresponding dispersion 
relations. The roots of the dispersion relations are chosen such that Re [q m ,c] > 0, thereby ensuring that 
the SPP wave decays in directions normal to the interface z = 0. 

For our numerical studies we chose e m = — 56 + 21i (which is the value for bulk aluminium at a free-space 
wavelength of 633 nm), in keeping with previous studies [3]. In Fig. 2, the real and imaginary parts of a are 
plotted against tj/L € (0, 1) and k$L € (0, 0.2). The real part of a increases uniformly as both the correlation 
length and the size parameter increase. Accordingly, a significant decrease in the phase speed of the SPP 
wave may be inferred as both L increases and rj increases. As the imaginary part of a similarly increases in a 
uniform fashion as both the correlation length and the size parameter increase, we deduce that attenuation 
of the SPP wave is significantly increased upon increasing both L and rj. 

Next we turn to the penetration of the SPP wave into the CTF in the +z direction. An inverse measure of 
this is provided by Re [q c ], which is plotted against r//L and k$L in Fig. 3. Thus, from Fig. 3 we deduce that 
SPP wave penetration into the CTF decreases as the correlation length increases and as the size parameter 
increases. 



4 Modified Kretschmann configuration 

The Kretschmann configuration [27] represents a realistic scenario for SPP wave excitation. A modification 
of this configuration is considered here [28] , guided by experimental views for launching surface waves [29] , 
wherein a fluid-infiltrated CTF of finite thickness occupies the region L m < z < and a thin metal film 
occupies the region < z < L m . The region z > Ly, is filled by a dielectric material with relative permittivity 
eg = n 2 , while the region z < is filled by a dielectric material with relative permittivity e<i- 

Suppose that a p-polarized plane wave propagates in the region z < towards the metal-coated CTF 
at an angle 6i nc £ (0,7r/2) to the +z axis. This incident plane wave gives rise to a reflected plane wave in 
the z < region and a transmitted plane wave in the z > Ls region. Thus, the corresponding electric field 
phasors in the regions above and below the metal-coated CTF may be expressed as [3, 26] 

!(— u x cos6» mc + u z sm9 inc ) exp [i (nx + ^Jl~dZ cos6» inc )] 
+r p (u x cos 8 lnc + u z sin 9 mc ) exp [i (nx - ^/e~dZ cos6 inc )] , z<0 ; (8) 

t p (—u x cos9tr + u z sin# tr ) exp [i (kx + sJTi (z — Ls) cos#tr)l i z > Ly, 
wherein the angles of incidence and transmission are related by 

y/7d sin 9 inc — ^/u sin 9 tr = a (9) 

and k = fcocr. By solving the associated boundary value-problem, the complex-valued reflection and trans- 
mission coefficients, namely r p and r t , are determined. 

We may identify the excitation of a SPP wave at the metal-CTF interface via the reflectance | r p | 2 or 
the absorbance 

A p = l-(\r p \ 2 +[3\t p \ 2 ), (10) 

wherein the scalar parameter 

^_ y/Ti Re [cos6> tr ] 

yJTd COS 9 inc 
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A sharp trough in graph of | r p | 2 (or sharp peak in the graph of A p ) versus 9i nc is the characteristic signature 
of SPP-wave excitation, provided that 9i nc is greater than the critical angle for total reflection in the absence 
of the metal film. In order to ensure the existence of this critical angle, in our numerical studies we chose 
Ed = 6.76, which is in the range of rutile. 

A representative plot of the reflectance | r p | 2 versus 9i nc is presented in Fig. 4 for k$L = 0.2 and n/L = 1. 
We set the thicknesses L m = 10 nm and = L m + 1000 nm. Also plotted in this figure is /3|i p | 2 . In the 
regime Q- mc > 36° (which corresponds to the total reflection regime in the absence of the metal film), the 
value of the absorbance A p can be inferred from that of the reflectance | r p | 2 since /3|t p | 2 is null- valued here. 
The rightmost local minimum in the graph of | r p | 2 (at 9i nc = 70.8°) indicates the excitation of a SPP wave. 
The other local minimums in the graph of | r p | 2 , which occur for 36° < 9i nc < 70.8°, represent waveguide 
modes [30] . Crucially, these waveguide modes are excited at different values 9t nc as the thickness of the CTF 
is varied, whereas the value of 9i nc which excites the SPP wave is independent of the CTF's thickness (as 
we confirmed in numerical calculations not presented here). 

Let us now focus on 9\ nc which denotes the value of 9i nc corresponding to the rightmost local minimum in 
the graph of | r p | 2 versus 8i nc ; i.e., it corresponds to the value of 9i nc at which the SPP wave is excited. The 
sensitivity of this angle 9\ nc to the correlation length and size parameter is demonstrated in Fig. 5, where 
®\nc i s plotted against n/L G (0, 1) and k n L £ (0,0.2). We see that 6\ nc increases uniformly as L increases 
and as 77 increases. In particular, the value of 9\ nc for L = 0.2/fco and 77 = L is approximately one degree 
greater than it is in the limits L — > and r\ — > 0. 

From the point of view of sensing applications, the shape of the local minimum at 9 mc = 9\ nc is an 
important consideration. To explore the effects of the correlation length and size parameter on the shape of 

the SPP trough, the second derivative d 2 (V P | 2 ^ /d9f nc is plotted against n/L G (0, 1) and k L € (0,0.2) in 

Fig. 6. We observe that the SPP trough becomes uniformly sharper as L increases, but changes relatively 
little as 77 increases. In particular, the SPP trough for L = 0.2/ kg and rj = L is approximately 2% sharper 
than it is in the limits L — > and r; — > 0. 

Lastly, we turn to the sensitivity of the SPP trough in the graph of |r p | 2 versus 9i nc to changes in the 
refractive index rig of the fluid which infiltrates the CTF. A convenient measure of this sensitivity is provided 
by the figure of merit (in degree/RIU 2 ) [3] 

n£2 - Tin 

where 9\ nc is regarded as function of rig. For ri£2 — 1-5 and ng\ = 1.0 the figure of merit is plotted against 
rj/L £ (0,1) and kgL £ (0,0.2) in Fig. 7. The sensitivity, as gauged by p, decreases uniformly as both 
the correlation length increases and as the size parameter increases. In particular, the sensitivity drops by 
approximately 4 degree/RIU as L and rj both increase from zero to 0.2/fc . 

5 Closing remarks 

The potential that infiltrated CTFs offer for SPP-based optical sensing has been demonstrated in previous 
theoretical studies using the Bruggeman homogenization formalism. We have further elucidated this matter 
by implementing a more sophisticated homogenization approach based on the extended second-order SPFT 
which takes into account the statistical distribution of the component materials and the size of the component 
materials. Specifically, our numerical investigations have revealed that as the correlation length increases 
and as the size parameter increases: 

• the phase speed of the SPP wave decreases and the SPP wave's attenuation increases; 

• the SPP wave's penetration into the CTF decreases; 

2 RIU = refractivc-indcx unit 
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• the angle of incidence required to excite the SPP wave in a modified Kretschmann configuration 
increases; 

• the sharpness of the SPP trough in the graph of reflectance versus angle of incidence increases; 

• and the sensitivity to changes in refractive index of the infiltrating fluid decreases. 

As well as directly shedding further light on SPP-based CTF sensors, these results also highlight the 
importance of using higher-order homogenization methods when considering sculptured thin films generally 
as platforms for optical sensing. For example, in principle a similar higher-order homogenization approach 
could be usefully employed for chiral sculptured thin films, which are also promising candidates for SPP-based 
optical sensing [7]. 
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Figure 1: The relative permittivity parameters e a 2, £fc2 and e C 2 of a fluid-infiltrated CTF, as computed using 
the extended second-order SPFT, plotted against rj/L € (0, 1) and koL G (0,0.2). The refractive index of 
the fluid infiltrating the CTF is fixed at ng = 1.5. „ 
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Figure 3: The real part of q c plotted against rj/L <G (0, 1) and k L S (0,0.2). 
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6» inc (degree) 

Figure 4: The reflectance |r p | 2 (red, solid curve) plotted against 6i nc (in degree), for the modified 
Krctschmann configuration with = 6.76, e m = —56 + 21i, eg = 2.25, L m = 10 nm, Ls = L m + 1000 
nm, r~i/L = 1.0 and koL = 0.2. Also plotted is the quantity j3 \t p \ 2 (blue, dashed curve). 




Figure 5: The angle 9\ nc (in degree) for the scenario represented in Fig. 4, plotted versus ?]/L e (0, 1) and 
k a L e (0,0.2). 
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Figure 7: The figure of merit p (in degree/RIU) plotted against n/L £ (0, 1) and koL £ (0,0.2) for = 1.5 
and nn = 1-0. 



12 



